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ABSTRACT 

We present a study of the resolved star-forming properties of a sample of distant 
massive (M > 10 n M o ) galaxies in the GOODS NICMOS Survey (GNS), based on 
deep Hubble Space Telescope imaging from the GOODS North and South fields. We 
derive dust corrected UV star formation rates (SFRs) as a function of radius for 
O j 45 massive galaxies within the redshift range 1.5 < z < 3 in order to measure the 

spatial location of ongoing star formation in massive galaxies. We find that the star 
formation rates present in different regions of a galaxy reflect the already existent 
stellar mass density, i.e. high density regions have higher star formation rates than 
lower density regions, on average. This observed star formation is extrapolated in 
several ways to the present day, and we measure the amount of new stellar mass that 
is created in individual portions of each galaxy to determine how the stellar mass 
added via star formation changes the observed stellar mass profile, the Sersic index 
and effective radius over time. We find that these massive galaxies fall into three 
broad classifications of star formation distribution: 1) Total stellar mass added via 
star formation is insignificant compared to the stellar mass that is already in place 
at high redshift. 2) Stellar mass added via star formation is only significant in the 
outer regions (R > lkpc) of the galaxy. 3) Stellar mass added via star formation 
is significant in both the inner (R < lkpc) and outer regions of the galaxy. These 
different star formation distributions increase the effective radii over time, which are 
on average a factor of ~ 16 ± 5% larger, with little change in the Sersic index (average 
An = — 0.9±0.9) after evolution. We also implement a range of simple stellar migration 
models into the simulated evolutionary path of these galaxies in order to gauge its effect 
on the properties of our sample. This yields a larger increase in the evolved effective 
radii than the pure static star formation model, with a maximum average increase 
of AR e ~ 54 ± 19%, but with little change in the Sersic index, An ~ —1.1 ± 1.3. 
These results are not in agreement with the observed change in the effective radius 
and Sersic index between z ~ 2.5 and z ~ obtained via various observational studies. 
We conclude that star formation and stellar migration alone cannot account for the 
observed change in structural parameters for this galaxy population, implying that 
other mechanisms must additionally be at work to produce the evolution, such as 
merging. 

Key words: 



1 INTRODUCTION 

One of the least understood aspects of galaxy evolution is 
the star formation rates in galaxies, how these vary across in- 
dividual galaxies, and influence galaxy properties. A key way 
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to address galaxy evolution directly is to understand how the 
nearby galaxy population was put into place and evolved 
from higher redshift galaxies, which we can now observe in 
near complete mass-selected samples up to 2 = 3 (e.g. Daddi 
et al. 2007; Conselice et al. 2011). One major finding of high 
redshift studies is that massive galaxies (M* > 1O 11 M0) 
have significantly smaller effective radii than low redshift 
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galaxies of similar mass (e.g. Daddi et al 2005; Trujillo et al. 
2006a, 2006b, 2007, 2011; Buitrago et al. 2008, 2011; Cimatti 
et al. 2008; van Dokkum et al. 2008, 2010; Franx et al. 2008 
; van der Wei et al. 2008; Damjanov et al. 2009; Carrasco et 
al. 2010; Newman et al. 2010; Szomoru et al. 2011; Weinzirl 
et al. 2011). 

Several physical processes have been proposed to ex- 
plain this strong size evolution within the massive galaxy 
population at z < 2. These can be divided into two distinct 
categories, external processes such as gas poor (dry) merg- 
ers (e.g. Khochfar & Silk 2006; Naab et al. 2009) and cold 
gas flows along cosmic web filaments (e.g. Dekel et al. 2009; 
Conselice et al. 2012 submitted) as a means for puffing up 
the stellar components of these massive galaxies, or internal 
processes such as adiabatic expansion resulting from stellar 
mass loss and strong AGN-fuelled feedback (e.g. Fan et al. 
2008, 2010; Hopkins et al. 2010A; Bluck et al. 2011). One 
process that has not been looked at in detail is the internal 
star formation distribution present within massive galaxies 
at high redshift, and whether this can account for the ob- 
served structural evolution. This can now be examined due 
to high resolution data from the GOODS NICMOS Sur- 
vey taken with the ACS and NICMOS-3 instruments on the 
Hubble Space Telescope (Conselice at al. 2011). 

We know that galaxies evolve significantly in stellar 
mass from observational studies showing that half of the 
stellar mass of present day galaxies is already in place by 
z ~ 1 (e.g. Brinchmann & Ellis 2000; Drory et al. 2004; 
Bundy et al. 2006; Perez-Gonzalez et al. 2008; Mortlock et 
al. 2011). The most massive galaxies (M* > 10 n M s ) ap- 
pear on average to have red rest-frame colours which we 
expect to see for galaxies dominated by old stellar popula- 
tions (Saracco et al. 2005; Labbe et al. 2006; Conselice et al. 
2007; Griitzbauch et al. 2011). However, Bauer et al. (2011b) 
show that ~ 80% of these massive red galaxies likely har- 
bour dusty star formation. This star formation over cosmic 
time could contribute large amounts of stellar mass to mas- 
sive galaxies, and depending on where this mass is created 
could affect their observable structural properties as they 
evolve. 

In the merger scenario, estimates for the total number of 
major mergers experienced by massive galaxies since z — 3 
is N m = 1.7 ± 0.5 (Bluck et al. 2009). This would imply an 
average stellar mass increase of, at best, a factor of two due 
to major mergers. However over the same period of time the 
effective radius of massive galaxies has increased on average 
by a factor of three for disk-like galaxies, and a factor of five 
for spheroid-like galaxies, effectively building up stellar mass 
in the outer regions of galaxies (see e.g. Buitrago et al. 2008, 
2011; Trujillo et al. 2007; Carrasco et al. 2010; van Dokkum 
et al. 2010). This additional stellar mass could arise from 
star formation already present at high redshift within these 
outer regions. 

To date studies have only looked at the total star forma- 
tion rates of these galaxies as a whole (e.g. Perez-Gonzalez 
et al. 2008; Cava et al. 2010; van Dokkum et al. 2010; Bauer 
et al. 2011; Griizbauch et al. 2011; Viero et al. 2011; Hilton 
et al. 2012), but have not examined the locations of the star 
formation within these galaxies. Thus, we combine the ob- 
served stellar mass profiles with the observed star formation 
profiles of high redshift massive galaxies in order to measure 
the effect stellar mass added via star formation over ~ 10 



Gyr has on different spatial regions, and to the total stellar 
mass profile. We also ascertain whether this star formation 
can account for the observed size evolution. 

Along with size evolution within the massive galaxy 
population there is also a change in overall morphology. The 
present day universe is populated by massive galaxies with 
early-type morphologies (e.g. Baldry et al 2004, Conselice et 
al. 2006). At earlier epochs, z > 1.5, observational studies 
have found that the massive galaxy population is dominated 
by galaxies with late-type morphologies (e.g. Buitrago et al. 
2011; Cameron et al. 2011; van der Wei et al. 2011; Weinzirl 
et al. 2011). This morphological shift can be seen via a 
change in Sersic index from low values, n <~ 2.5 denoting a 
possible late-type morphology, to high values, n >~ 2.5 de- 
noting a possible early-type morphology. In the hierarchical 
model of galaxy evolution there are many methods that can 
drive morphological evolution. These methods include in situ 
star formation producing disk-like systems (e.g. Dekel et al. 
2009; Oser et al. 2010; Ricciardelli et al. 2010; Wuyts et al. 
2010; Bournaud et al. 2011), and/or mergers with satellite 
galaxies producing a more spheroid-like system (e.g. Khoch- 
far & Silk 2006; Hopkins et al. 2009; Feldmann et al. 2010; 
Oser et al. 2010). We therefore also investigate how in situ 
star formation over cosmic time changes the Sersic index of 
the massive galaxies, and ascertain whether this process can 
account for the observed morphological changes. 

This paper is set out as follows: Section 2 discusses the 
GOODS NICMOS Survey, the galaxy sample, and how the 
data used in this paper was obtained. Section 3.1 examines 
the stellar mass radial density distributions of the massive 
galaxies. Section 3.2 describes how the stellar mass density 
added via star formation is calculated. In Section 3.3 we 
examine the evolved galaxy profiles. Section 4.1 presents the 
findings of how the structure and size of the massive galaxies 
is altered by star formation. In Section 4.2 we introduce a 
simple stellar migration model to the stellar mass added 
by star formation in order to gauge the effect this has on 
structures and sizes. Section 5 and 6 contain the discussion 
and summary of our findings, respectively. Throughout this 
paper we assume Qm = 0.3, J1a = 0.7 and Ho — 70 km 
s" 1 Mpc" 1 . AB magnitudes and a Salpeter IMF are used 
throughout. 



2 DATA AND ANALYSIS 

In this section we describe the survey we use in this study, 
the GOODS NICMOS Survey (GNS), as well as the mea- 
surements of photometric redshifts, stellar masses, rest- 
frame colours and star formation rates for our galaxies. 

2.1 The GOODS NICMOS Survey 

The data used in this paper is obtained through the GOODS 
NICMOS Survey (GNS). The GNS is a 180 orbit Hubble 
Space Telescope survey consisting of 60 single pointings with 
the NICMOS-3 near-infrared camera, with an imaging depth 
of three orbits per pointing (Conselice et al. 2011). 

These pointings were optimised to contain the maxi- 
mum number of massive galaxies (M > 10 1 Mq) in the 
redshift range 1.7 < z < 3, identified in the two GOODS 
fields by their optical-to-infrared colours (see Conselice et al 
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2011). The survey covers a total area of about 45 arcmin 2 
with a pixel scale of ~ 0.1 arcsec/pixel, corresponding to 
~0.9 kpc at the redshift range of interest (1.5 < z < 3). The 
target selection, survey characteristics and data reduction 
are fully described in Conselice et al. (2011). Other analyses 
of the GNS data set can be found in Bauer et al. (2011), 
Bluck et al. (2010), Buitrago et al. (2008), Griitzbauch et 
al. (2011) and Mortlock et al. (2011). 

The GNS has a 5a limiting magnitude of Hab = 
26.8, which is significantly deeper than ground based near- 
infrared imaging of the GOODS fields carried out with e.g. 
ISAAC on the VLT, which reaches a 5a depth of Hab = 24.5 
(Retzlaff et al 2010). Sources were extracted from the NIC- 
MOS _Hi6o-band image and matched to the optical HST- 
ACS bands B,V,i and z, which are available down to a 
AB limiting magnitude of B = 28.2. The matching is done 
within a radius of 2 arcsec, however the average separation 
between optical and _Hi6o-band coordinates is much better 
with ~ 0.28 ±0.4 arcsec, roughly corresponding to the NIC- 
MOS resolution (see also Bauer et al. 2011b). 

The photometric catalogue covering the BvizH bands 
comprises 8298 galaxies, and is used to compute photometric 
redshifts, rest-frame colours and stellar masses described in 
the following sections (see also Conselice et al. 2011 for more 
details). Along with this, each galaxy has imaging data in 
the Bviz ACS bands (Giavalisco et al. 2004). 

Within our NICMOS fields we find a total of 81 galaxies 
with stellar masses larger than 10 11 Mq with photometric 
and spectroscopic redshifts in the range 1.5 < z < 3. This 
-Hi6o-band sample of massive galaxies is reduced to 52 due to 
optical band non-detections where we are unable to calculate 
accurate ultraviolet (UV) dust extinction corrections (Bauer 
et al. 2011). The sample is further reduced to 45 galaxies due 
to removing those galaxies with Sersic fits to the Hieo light 
profiles with high uncertainties (Buitrago et al. 2008). We 
examined the UV surface brightness profiles of the excluded 
galaxies in the bands they were detected in, and found them 
to be consistent with the profiles of the remaining galaxies. 
Figure [T] shows the zsso and i/i60 band images of the 45 
galaxies used in this study. 

2.2 Redshifts 

Where possible, we use spectroscopic redshifts published in 
the literature for our GNS galaxies, otherwise we use our 
own measured photometric redshifts. Spectroscopic redshifts 
of sources in the GOODS-N field were compiled by Barger et 
al. (2008) , whereas the GOOD-S field spectroscopic redshifts 
are taken from the FIREWORKS compilation (Wuyts et al. 
2008). In the full GNS sample, there are 537 spectroscopic 
redshifts for sources in GOODS-N and 369 in GOODS-S. 
In the massive galaxy sample used in this paper there are 
however only six galaxies with spectroscopic redshifts. 

Photometric redshifts are therefore crucial for this 
study. These photo-zs were obtained by fitting template 
spectra to the BvizH photometric data points using the HY- 
PERZ code (Bolzonella et al. 2000). The method is described 
in more detail in Griitzbauch et al. (2011b). The synthetic 
spectra used by HYPERZ are constructed with the Bruzual 
& Chariot evolutionary code (Bruzual & Chariot 1993) rep- 
resenting roughly the different morphological types of galax- 
ies found in the local universe. We use five template spectra 



corresponding to the spectral types of E, Sa, Sc and Im, 
as well as a single starburst scenario. The reddening law is 
taken from Calzetti et al. (2000). HYPERZ computes the 
most likely redshift solution in the parameter space of age, 
metallicity and reddening. The best fit redshift and corre- 
sponding probability are then output together with the best 
fit parameters of spectral type, age, metallicity, ^l^and sec- 
ondary solutions. 

To assess the reliability of our photometric redshifts 
we compare them to available spectroscopic redshifts in 
the GOODS fields. We matched the two catalogues to our 
photometric catalogue with a matching radius of 2 arcsec, 
obtaining 906 secure spectroscopic redshifts. The reliabil- 
ity of photometric redshifts measures we use is defined by 
Az/(1 + z) = (z apec - z photo )/(l + z apec ). In the follow- 
ing we compare the median offset from the one-to-one re- 
lationship between photometric and spectroscopic redshifts, 
(Az/(1 + z)), and the RMS scatter around this relation, 
v&z/ii+z) = 0.061. We then investigate the performance of 
HYPERZ at different redshifts, at low redshift (z < 1.5) 
and at 1.5 ^ z ^ 3, which is the redshift range of the galaxy 
sample we use. For the high redshift complete sample we 
obtain an average offset {Az/(1 + z)) — 0.06 and a RMS 
of (Taz/(i+ 2 ) = 0.1, with a fraction of catastrophic outliers 
of 20%, where catastrophic outliers are defined as galaxies 
with |Az/(l + z)| > 0.3, which corresponds to ~ 3 times the 
RMS scatter. This error has been folded into the results. 

2.3 Stellar Masses and e- folding Star formation 
Times 

Stellar masses and rest-frame colours of our sample are de- 
termined from multicolour stellar population fitting tech- 
niques using the same catalogue of five broad band data 
points used to determine photometric redshifts for all GNS 
galaxies. A detailed description of how stellar masses and 
rest-frame (U — B) colours are derived can be found in Con- 
selice et al. (2011) and Griitzbauch et al.(2011), and is sum- 
marised in the following. 

To calculate these masses and colours we construct a 
grid of model spectral energy distributions (SEDs) from 
Bruzual & Chariot (2003) stellar population synthesis mod- 
els, assuming a Salpeter initial mass function and a varying 
star formation history, age, metallicity and dust extinction. 
The star formation history is characterised by an exponen- 
tially declining model of the form 

SFR(t) = SFR x e~ t/T . (1) 

The parameters in Equation 1 are varied over a wide 
range of values within the ranges; r = 0.01 to 10 Gyr - with 
values (all in Gyr): 0.01 0.025 0.05 0.08 0.1 0.12 0.65 1.21 
1.37 1.73 1.99 2.52 2.71 2.84 2.99 3.09 3.34 3.7 4.23 4.33 4.45 
4.69 5.16 5.23 5.47 5.68 5.83 6.21 6.61 6.95 7.03 7.27 7.37 
7.95 8.38 8.76 8.80 8.94 9.57 9.80 9.88, and the time since the 
onset of star formation ranging from t = to 10 Gyr, with a 
condition that ages are not older than the universe itself at 
the redshift of observation. The dust content is parametrised 
by the V-band optical depth with values tv = 0.0, 0.5, 1,2 
and the metallicity ranges from 0.0001 to 0.05 (Bruzual & 
Chariot 2003). 

The magnitudes obtained from the model SEDs are fit 
to the observed photometric data of each galaxy using a 



© 2012 RAS, MNRAS 000.[Tlfl8l 



4 Ownsworth et al. 



Bayesian approach. A likelihood distribution stellar mass, 
age and absolute magnitude at each possible star formation 
history is computed for each galaxy. The stellar mass is de- 
termined based on this distribution, where the most likely 
stellar mass produces a peak in the distribution, and the 
uncertainty is the width. The final error, as a result of the 
models used, lie within the range of 0.2 to 0.3 dex. 

We use negative r models in this paper to fit our stel- 
lar masses, although other forms of the star formation are 
included later in the paper when considering the addition 
of stellar mass, including constant star formation rate mea- 
sures, and an exploration of possible forms of the star forma- 
tion, including those which are maximal and would exceed 
the observed stellar mass density evolution (§5.4.1). While 
there is some evidence that the star formation history ac- 
tually increases from z — 8 to z = 3 (e.g., Papovich et al. 
2011), there is also evidence that at redshifts lower than 
this, and particularly for high mass galaxies, that the star 
formation rate is starting to decline (e.g., Conselice et al. 
2007; 2011). We investigate this in detail for our sample by 
holding at a constant co-moving volume, as Papovich et al. 
(2011) does, and seeing how our star formation rate changes 
for the same co-moving density. Doing this, we find that the 
star formation declines over our epoch using this method, 
although the star formation history at redshifts z > 3 is 
more complicated than this. However, this does show that 
our values of r that we use here are mimicking the form of 
the empirical star formation history. 

It is possible that the stellar masses are an over esti- 
mate due to the poor treatment of the TP-AGB phase in 
a star's life. The effects of this phase are less important at 
the rest frame wavelengths used in this study, especially in 
the infra-red -Hi60 band. Using newer models by Bruzual & 
Chariot (2010) which have an improved treatment of the 
TP-AGB phase we find that this lowers the stellar masses of 
the massive galaxy sample by < 0.07 dex. This effect from 
the new models is smaller than the stellar mass error, and 
the effects of cosmic variance, and is therefore negligible. 
Table [T] contains the full list of values for all variables used 
in this study. 

2.4 Star Formation Rates 

The star formation rates (SFRs) used in this paper are mea- 
sured from rest-frame UV luminosities, using the methods 
described in Bauer et al. (2011). The rest-frame UV provides 
a direct measurement of ongoing SFR, since the UV lumi- 
nosity is directly related to the presence of young and short- 
lived stellar populations produced by recent star formation. 
However, UV light is very susceptible to dust extinction and 
a careful dust-correction has to be applied. The correction 
we use here is based on the rest-frame UV slope (Bauer at 
al. 2011). We briefly describe the method in the following. 

We determine the SFR uv ,uncorrected from the observed 
optical ACS zsso-band flux density (with a 5 a limit of 27.5 
in the AB system) spanning wavelengths of 2125 - 3400A 
for z — 1.5 - 3 galaxies. After applying an SED based k- 
-correction using the IDL KCORRECT package (Blanton & 
Roweis 2007, v4.2). This is done by using the full SEDs of 
these galaxies. This result of this is a k-correction at rest- 
frame UV wavelengths ~ 2800A which we use throughout 
this paper. 



To measure the SFR we first derive the UV luminosity 
of the massive galaxies, then use the Kennicutt (1998) con- 
version from 2800A luminosity to SFR assuming a Salpeter 
IMF: 

SFRuviMoyr- 1 ) = 1.4 x 10- 28 L 2 80o(ergss- 1 Hz" 1 ) (2) 

Before dust extinction is taken into account we find a limit- 
ing SFRuv,obs = O.3±O.lM0yr _1 at z = 1.5, and a limiting 
SFR UVt0bs = 1.0 ± O.3M yr _1 at z = 3. The errors quoted 
here take into account photometric errors and the conver- 
sion from a luminosity. The error for individual SFRs are 
around 30%. This error is dominated by the dust correction 
applied. 

We compare our total integrated SFR for each galaxy in 
this sample to the same sample used in Bauer et al. (2011) 
which included both SED determined and UV SFR. We find 
on average that our total SFRs are slightly higher due to the 
use of a larger aperture but within the quoted error. 

Several studies (e.g. Bauer et al. 2011; Reddy et al. 
2012) have found that in when comparing the IR derived 
SFRs plus UV derived SFRs (SFRir+uv) against dust cor- 
rected UV SFRs (SFRuv,corr) that the SFR IR+UV is on 
average a factor of 3 larger than SFRuv,corr- This overesti- 
mation has been seen in other studies looking at luminous 
galaxies (e.g. Papovich et al. 2007). Results from the Her- 
schel Space Telescope (e.g. Elbaz et al 2010; Nordon et al 
2010; Hilton et al. 2012) suggest that at z > 1.5, the 24 
/jm flux may overestimate the true SFR due to a rise in 
the strength of polycyclic aromatic hydrocarbon (PAH) fea- 
tures, changes in the SEDs, or AGN contamination. 

Recent work on the same sample of massive galaxies 
used in this paper by Hilton et al. (2012) using Herschel 
Space Telescope data and new fitting methods of the same 
sample of massive galaxies found that the scatter of these 
results can be reduced to < la between the IR+UV and UV 
corrected SFRs by taking into account new templates in the 
FIR that account for these issues. We however are unfortu- 
nately forced to only use the SFRuv,corr in this paper since 
the Spitzer IR 24^m and Herschel images are not resolved. 

2.5 Dust Corrections 

To obtain reliable star formation rates in the rest-frame ul- 
traviolet, we need to account for the obscuration due to dust 
along the line of sight. Meurer et al. (1999) found a corre- 
lation between attenuation due to dust and the rest-frame 
UV slope, P, for a sample of local starburst galaxies (where 
F\ ~ A 8 ). More recent studies of local galaxies using the 
Galaxy Evolution Explorer (GALEX) near-ultraviolet band 
show that the UV slope from the local starburst relation 
can be used to recover the dust attenuation of moderately 
luminous galaxies at z ~ 2 (Buat et al. 2005; Seibert et al. 
2005; Reddy et al. 2010). 

A method for determining dust extinction uses the red- 
dening parameter extracted from the best-fitting SED tem- 
plate as described in §2.3. We fit sets of template stellar 
population synthesis models to derive the stellar masses 
(Griitzbauch et al. 2011; Conselice et al. 2011). This method 
has some limitations when using it to correct for dust as this 
approach assumes that the UV slope is due to dust reddening 
instead of other sources, such as evolved stellar populations 
(Griitzbauch et al. 2011). 
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We apply a method for determining a UV dust attenu- 
ation, A2800, in terms of the UV slope. The UV slope is de- 
termined using an SED-fitting procedure described in Bauer 
et al. (2011). To summarise, we fit an SED to the multi- 
wavelength observations from optical-to-infrared. The SEDs 
obtained for all sources in the GOODS fields were fit with 
stellar population synthesis models. The best fitting tem- 
plates were then used to obtain a synthetic estimate of the 
UV emission at 1600 A and 2800 A. From the model-derived 
UV luminosities at 1600A and 2800A we calculate the spec- 
tral slope, /3. The Calzetti et al. (2000) law is then used to 
derive A2800 from the UV spectra slope, which we apply to 
the UV-derived star formation rates. Using this method we 
find an average extinction value of ^2800 = 3.2 ± 1.0 magni- 
tudes for our sample. 

Bauer et al. (2011) find in a comparison between an 
SED determined, and an observationally derived A2800, that 
the values obtained from these two methods are in relatively 
good agreement for M* > lO n M0 across the whole redshift 
range with an average offset of SA2S00 = 0.86. This corre- 
sponds to a ~ 27% error in the average dust attenuation 
and this is folded into the following results. 



3 STELLAR MASS DENSITY PROFILES 
3.1 Stellar Mass Radial Density Distributions 

We construct our sample galaxy's stellar mass density 
distributions by examining the distribution of the -ffi6o, rest- 
frame optical light profiles for our sample (for 1.5 < z < 3 
this corresponds to rest A = 640 — 400nm). We base this 
determination on the Sersic fits to the light distribution. 
A detailed description of how the Sersic indices were mea- 
sured can be found in Buitrago et al. (2008, 2011) and is 
summarised in the following. 

The Sersic profiles were measured using GALFIT (Peng 
et al. 2002, 2010). GALFIT uses r 1/n 2D models of the form 
(Sersic 1968): 

E(r) = S e x exp(-b n [(R/R c ) 1/n - 1]) (3) 

Where R e is the effective radius of the galaxy, S e is the 
surface brightness at R e , n is the Sersic index and b n = 
1.9992n - 0.3271. This model is convolved with the Point 
Spread Function (PSF) of the images, and GALFIT deter- 
mines the best fit by comparing the convolved model with 
the observed galaxy surface brightness distribution using 
a Levenberg-Marquardt algorithm to minimise the \ 2 °f 
the fit. We used single Sersic models to compare our size 
estimations with previous work. Neighbouring galaxies are 
masked out before the fitting, and in the case of overlapping 
isophotes the objects are fit simultaneously. Due to varia- 
tions of the shape of the NICMOS-3 PSF in our images, we 
select five non- saturated bright stars to sample the PSF 
within our imaging and with which to gauge the accuracy 
of the parameter measurements. The structural parameters 
of each galaxy are measured five times for each unique star. 
The uncertainty (la) on the structural parameters due to 
changes in the PSF is ~ 15% for the effective radius r e , and 
~ 20% for the Sersic index n. 

A concern when measuring sizes and Sersic indices at 
high redshift is surface brightness dimming which in prin- 



cipal could bias our measured sizes. Previous studies have 
examined this issue, and have conducted many simulations 
in order to check the importance of surface brightness dim- 
ming in HST observations (e.g. Trujillo et al. 2006a, 2007; 
Buitrago et al. 2011). In Appendix A of Buitrago et al. 
(2011) one can find the descriptions of the extensive sim- 
ulations we conducted in order to asses the reliability of the 
galaxy structural parameters within GNS. The median ob- 
servable characteristics of our massive galaxies (Hab = 22.5, 
n~2, r e ~2 kpc) allow us to retrieve their structural proper- 
ties without any significant bias. However it is worth noting 
that, for individual galaxies, the parameters are not as well 
constrained for galaxies which display higher Sersic indices. 

Using the total stellar mass for the galaxies and the 
H-band light (see §2.3) we equate the total stellar mass to 
the total rest-frame optical light received from the individ- 
ual galaxies. We then convert the _ffi6o-band Sersic profile 
to a stellar mass profile. With the total initial stellar mass 
defined by: 

M, = p e / exp(-b n [(R/R c ) 1/n - l])27rRdR (4) 
Jo 

The radius within which the total stellar mass is contained, 
Rmax, is taken to be 20kpc in all cases. The effective radius, 
R e and Sersic index, n, come from the Hi6o-band Sersic 
profile as described in this section. From this, the stellar 
mass density at the effective radius p e is calculated, and the 
full stellar mass density profile is constructed via: 

Pobserved{R) = pe X exp(-b n [(R/Re ) 1 '° - 1]) (5) 

with the implicit assumption that the mass to light ratio is 
constant over the galaxy, as used in other works studying 
surface brightness profiles (e.g. Szomoru et al. 2011). 

3.2 Stellar Mass Density Added Via Star 
Formation 

We measure star formation profiles for our sample using 
the IRAF program ellipse by fitting a series of isophotal el- 
lipses to the 2850-band data with the Hi6o-band determined 
centre of the massive galaxies. This isophotal fitting returns 
the 2850-band flux binned in a series of increasing radii. This 
is then converted to a dust corrected star-formation rate in 
each radius bin via the procedure described in §2.5. Exam- 
ples of such surface brightness profiles are shown in Figure 

m 

The total galaxy magnitudes obtained from this isopho- 
tal fitting are checked against a previous catalogue of the 
^850-band magnitudes for these galaxies in Bauer et al (2011) 
measured with SExtractor, and are found to be consistent. 

In order to measure how this star formation affects the 
local mass density of the galaxy we simulate the amount 
of stellar mass added via star formation in each radius bin 
by assuming that the same global star formation history we 
use in §2.3 applies throughout. We also apply several other 
star formation histories that these galaxies could experience 
such as constant SFR to z = 0, constant SFR to z — 1.5 
and variations on the derived tau model. These are disused 
in §5.4.2. We find very similar results as discussed below for 
the tau models. 

The exponentially declining model of star formation 
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Figure 1. Images in the 2850 band (left) and Hieo band (right) of the 45 massive galaxies in this study. The ID number of each galaxy 
is shown in the lower right hand corner of the zg$Q band image with the corresponding image of the same region in the Hiqq band to the 
right. All images are 2.5 by 2.5 arcsecond cutouts for galaxies between 1.5 < z < 3 centred on the -H160 band detection. The properties 
of each galaxy are listed in Table 1 . 



uses the observed dust-corrected rest-frame UV star forma- 
tion as the initial rate, SFRq, and the values of the e- folding 
time, r, which we obtained from the M„ fitting (see §2.3). 
To obtain the total amount of stellar mass added via star 
formation, Msf, we integrate Equation 1 over time from the 
total look back time derived from the redshift of the galaxy, 
ranging from ~ 9.7 Gyr at z — 1.7 to ~ 11.5 Gyr at z = 3 
to the present day. We experimented with evolving the mas- 
sive galaxies only until z — 1 but found very similar results 
as the evolution to z = as the majority of the evolution in 
both size and structural properties of these massive galaxies 
seems to occur within the first ~ 2 Gyr of our simulation. 

In Figure [3] we show the change in the total stellar 
mass of each of the galaxies within our sample as measured 
through the SFR. We find that the total stellar mass on 
average increases by 90% via this modelled star formation. 
The evolved total stellar masses of our galaxies do not ex- 



ceed constraints placed upon the observed total stellar mass 
evolution from other studies (e.g. Cole et al. 2001; Bell et al. 
2003; Conselice et al. 2007; Brammer et al. 2011; Mortlock et 
al, 2011). Brammer et al. (2011) show that the total stellar 
mass growth for massive galaxies from z ~ 2 to is on the 
order of 100%. We see that there is an anticorrelation be- 
tween the original stellar mass and the evolved stellar mass, 
with some of the lower mass galaxies increasing substan- 
tially in stellar mass, while the higher mass galaxies have a 
much smaller change in mass over cosmic time. This is a sign 
of galaxy downsizing (e.g. Cowie et al. 1996; Bundy et al. 
2006), such that the most massive galaxies are less affected 
by star formation at z < 3 than the lower mass galaxies. 

This total stellar mass added via star formation is 
converted to a stellar mass projected density via, psn = 
Msf I 'A an , were A an is the area of the annulus the star for- 
mation is contained within. From this we construct a new 
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ID: 3629 




R (Kpc) 

(a) 



R (Kpc) 

(b) 



Figure 2. Example rest frame UV surface brightness profiles from ACS zsso-band imaging, (a) Galaxy ID: 999, Initial stellar mass: 
1.5 X 10 Mq, Rest frame optical effective radius: 2.0kpc, Rest frame optical Sersic index:n = 1.42, SF growth classification: Inner 
SF growth, (b) Galaxy ID: 3629, Stellar mass: 2.6 X IO^Mq, Rest frame optical effective radius: 1.8kpc, Rest frame optical Sersic 
index:n = 1.26, SF growth classification: Non-significant SF growth. 



stellar mass profile by including the stellar mass added via 
star formation to the original profile via: 

M,(R,t) = M«(R,t = 0) + SFR (R) [ e~ t/T dt (6) 

Jo 

Where M* (R, t = 0) is the initial stellar mass at radius R, 
SFRo(R) is the observed initial SFR at the same radius 
and M* (R, t) is the stellar mass at radius R after time t. We 
consider other cases of SF histories in §5.4. 
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Figure 3. Total stellar mass before (Original M*) and after evo- 
lution (Evolved M,) from the derived tau model of star formation 
evolution. The black circle represent the non-significant star for- 
mation growth galaxies. The blue squares represent the inner star 
formation growth galaxies. The red triangles represent the outer 
star formation growth galaxies (see §3.3). 



3.3 Profiles 

In this section we examine the profiles of the stellar mass 
density distributions at high redshift, the stellar mass den- 



sity added via star formation, and for the combination of the 
two - an evolved stellar mass density profile for each of the 
45 massive galaxies in our sample. The evolved profiles are 
then fit with a new Sersic profile. The Sersic profiles for the 
evolved stellar mass profiles are obtained by the best fitting 
Sersic (1968) function to the new profile, 

p(R) = p e x exp(-b n [(R/R c ) 1/n - 1]) (7) 

We find that the galaxies in our sample can be classified into 
three distinct groups based on the location of star forma- 
tion regions and the effect they have on our sample galaxy's 
evolved stellar mass density profile. 

To examine the results we first divide the galaxies into 
two regions. An inner region, at R ^ lkpc with an observed 
initial stellar mass density in the inner region, p bserved,inner 
and a stellar mass density added via star formation in the 
inner region, psF,inner- An outer region, R > lkpc with a 
observed stellar mass density, pobserved, outer and a stellar 
mass density added via star formation in the outer region, 
Psf, outer- We chose lkpc as the boundary for our inner re- 
gion based on stellar surface brightness comparisons at high 
and low redshifts (e.g. Hopkins et al. 2009; Carrasco et al. 
2010; Szomoru et al. 2011). We discuss the three types be- 
low. 

Non-significant Star Formation Growth (NG) : This 
category is for galaxies in which the stellar mass density 
added via star formation is smaller then the galaxy's initial 
stellar mass density present over both the inner and outer 

regions. Pobserved,inner ^ PS F,inner and Pobserved, outer ^* 
PSF, outer- 

Outer Star Formation Growth (OG) : In this category 
the stellar mass density added via star formation is greater 
than the initial stellar mass density present in the outer 
region, but the initial stellar mass density in the inner re- 
gion is greater than the stellar mass density added via star 

formation, pobserved,inner ^ PS F,inner but pobserved, outer *C 
PSF, outer 
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ID 


SF Class 


SFR 




Mass 




R e 




1 ) 




(yr) 


43 


OG 


126.7±34.2 


1.79 


11 


.0±0.4 


3.8±1.0 


1.8±Q. 


3 


2.5±0.8 


6.5 


X 


10 8 


77 


NG 


173.8±46.9 


2.33 


11 


1±0.2 


3.5±0.9 


3.2±0. 


1 


1.4±0.1 


1.2 


X 


ur 


158 


NG 


28.8±7.8 


1.84 


11 


.2±0.4 


2.0±0.5 


1.4±0. 


3 


2.0±2.0 


6.5 


X 


10 8 


227 


NG 


23.1±6.2 


2.48 


11 


.2±0.3 


0.8±0.2 


2.1±Q. 


1 


2.1±0.1 


1.2 


X 


10 s 


840 


NG 


87.8±23.7 


2.31 


11 


.1±0.3 


2.4±0.6 


1.6±Q. 


2 


2.3±0.5 


1.2 


X 


10 8 


856 


NG 


257.8±69.6 


2.32 


11 


.2±0.2 


2.6±0.7 


1.7±0. 


1 


3.7±0.3 


1.2 


X 


10 s 


860 


OG 


367.0±81.0 


1.79 


11 


.2±0.3 


3.9±1.1 


3.8±0. 


3 


3.5±0.6 


6.5 


X 


10 a 


999 


IG 


672.6±181.6 


1.58 


11 


.2±0.4 


4.Q±1.1 


2.0±0. 


1 


1.4±0.3 


1.4 


X 


10 9 


1129 


NG 


464.7±125.5 


2.61 


11 


.3±0.4 


5. Oil. 3 


3.1±0. 


1 


0.2±0.1 


1.2 


X 


10 s 


1394 


NG 


208.4±56.3 


2.29 


11 


.5±0.3 


3.5±0.9 


3.3±0. 


1 


l.lrtO.l 


1.2 


X 


10 s 


1533 


NG 


97.4±26.3 


2.45 


11 


.6±0.3 


1.8±0.5 


1.3±Q. 


1 


0.9±0.5 


6.5 


X 


Kr 


1666 


NG 


244.5±66.0 


*1.76 


11 


.9±0.4 


4.7±1.3 


6.2±Q. 


1 


0.6±0.1 


1.2 


X 


10 8 


1768 


NG 


314.7±85.0 


2.22 


11 


.4±0.3 


3.6±1.0 


1.4±0. 


1 


2.1±0.6 


1.2 


X 


10 <s 


1888 


NG 


133.4±36.0 


2.75 


11 


.2±0.4 


3.5±0.9 


2.0±0. 


5 


5.0±1.2 


1.2 


X 


10* 


2083 


NG 


160.6±43.4 


2.31 


11 


.1±0.4 


3.5±0.9 


1.2±Q. 


2 


1.3±0.6 


1.2 


X 


Kr 


2411 


OG 


241.7±65.3 


2.09 


11 


.0±0.3 


3.8±1.0 


5.2±Q. 


1 


0.6±0.1 


6.5 


X 


10 8 


2564 


NG 


127.8±34.5 


2.10 


11 


.1±0.3 


3.5±0.9 


1.4±0. 


1 


1.8±0.3 


1.2 


X 


10 s 


2667 


OG 


284.8±76.9 


1.79 


11 


.2±0.4 


3.5±0.9 


2.0±0. 


1 


3.5±0.9 


6.5 


X 


10 s 


2678 


NG 


51.8±14.0 


2.30 


11 


.1±0.3 


1.2±0.3 


1.4±Q. 


1 


0.8±0.1 


1.2 


X 


U.r 


2798 


NG 


196.2±53.0 


1.72 


11 


.6±0.3 


3.5±0.9 


4.5±0. 


6 


4.6±0.2 


6.5 


X 


l() lS 


3629 


NG 


140.6±37.0 


2.17 


11 


.4±0.2 


2.5±0.7 


1.8±0. 


2 


1.3±0.2 


6.5 


X 


10 8 


3766 


NG 


76.0±20.5 


1.87 


11 


.2±0.2 


2.2±0.6 


2.8±0. 


1 


1.9±0.3 


1.4 


X 


10 9 


3818 


NG 


255.9±69.1 


1.82 


11 


.5±0.4 


3.5±0.9 


3.7±Q. 


4 


5.1±1.3 


6.5 


X 


1() 8 


3822 


NG 


115.3±31.1 


2.41 


11 


.1±0.5 


3.5±0.9 


1.9±0. 


1 


2.1±0.5 


1.2 


X 


If) 8 


4033 


OG 


81.5±22.0 


1.72 


11 


.2±0.3 


1.0±0.3 


6.1±0. 


8 


1.2±0.1 


1.7 


X 


10" 


4036 


OG 


134.3±36.3 


1.72 


11 


.0±0.3 


3.5±0.9 


1.7±0. 


1 


7.5±2.0 


6.5 


X 


■i nR 

10 s 


4315 


NG 


231.1±62.4 


2.85 


11 


.0±0.4 


4.Q±1.1 


0.9±0. 


5 


1.8±1.3 


1.2 


X 


It) 8 


4471 


OG 


401.0±108.3 


*2.29 


11 


.2±0.4 


3.5±0.9 


2.8±0. 


1 


1.4±0.1 


6.5 


X 


i ™R 
It) 8 


4557 


NG 


232.6±62.8 


2.09 


11 


.3±0.4 


3.1±0.8 


4.6±0. 


1 


1.4±0.1 


1.2 


X 


10 8 


4706 


NG 


236.2±63.8 


*2.35 


11 


.1±0.2 


3.4±0.9 


0.9±0. 


1 


1.7±1.2 


1.2 


X 


10 s 


4737 


NG 


34.4±9.2 


2.52 


11 


.0±0.3 


1.3±0.4 


1.1±Q. 


3 


3.3±3.7 


1.2 


X 


t /iR 
It) 8 


4882 


IG 


448.4±121.1 


1.67 


11 


.1±0.3 


4.2±1.1 


1.2±Q. 


1 


2.2±0.2 


6.5 


X 


i /iR 
It) 8 


5282 


IG 


276.8±74.7 


1.64 


11 


.0±0.3 


3.5±0.9 


4.9±0. 


1 


1.1±0.1 


2.7 


X 


10 9 


5764 


NG 


170.2±41.6 


2.54 


11 


.5±0.3 


3.5±0.9 


1.7±0. 


1 


2.9±0.4 


1.2 


X 


■i nfi 

10 8 


6035 


NG 


208.6±56.3 


1.60 


11 


.3±0.4 


3.5±0.9 


2.0±0. 


1 


4.2±0.6 


6.5 


X 


If) 8 


6114 


NG 


108.0±29.2 


1.96 


11 


.3±0.3 


2.5±0.7 


1.5±Q. 


6 


0.3±0.3 


6.5 


X 


10 s 


6220 


IG 


558.3±150.7 


1.71 


11 


.0±0.4 


5.1±1.4 


3.8±0. 


1 


1.5±0.1 


6.5 


X 


If) 8 


6267 


OG 


136.2±36.8 


*1.54 


11 


.4±0.2 


1.6±0.4 


2.3±0. 


1 


2.3±0.5 


1.7 


X 


10 () 


6287 


NG 


102.9±27.8 


1.84 


11 


.0±0.4 


3.5±0.9 


2.5±0. 


1 


2.4±0.4 


6.5 


X 


If) 8 


6514 


NG 


199.6±53.9 


2.49 


11 


.2±0.3 


3.2±0.9 


2.7±0. 


2 


2.5±0.1 


1.2 


X 


10 s 


6584 


NG 


180.4±48.7 


1.62 


11 


.2±0.3 


3.5±0.9 


6.2±0. 


1 


1.1±0.1 


1.2 


X 


10 8 


7072 


OG 


101.0±27.3 


1.74 


11 


.1±0.2 


1.9±0.5 


1.4±0. 


1 


2.3±0.4 


1.4 


X 


10 9 


7475 


OG 


377.9±102.0 


*1.61 


11 


.2±0.3 


3.5±0.9 


4.9±0. 


2 


2.1±0.2 


6.5 


X 


10 8 


8140 


OG 


319.6±86.3 


*1.90 


11 


.4±0.4 


3.0±0.8 


1.8±Q. 


1 


2.7±0.5 


6.5 


X 


If) 8 


8214 


OG 


817.0±220.6 


2.05 


11 


.3±0.2 


3.8±1.0 


1.6±Q. 


1 


1.7±0.2 


6.5 


X 


10 8 



Table 1. (col. 1) ID number of the galaxy; (col. 2) The classification of the galaxy based on the location of the star formation (see §3.3); 
Non-significant star formation growth (NG), outer star formation growth (OG) and inner star formation growth (IG) ; (col. 3) Total 
observed UV star formation rate in solar masses per year ; (col. 4) Best redshift of the object, spectroscopic redshifts are denoted by * 
; (col. 5) Stellar Mass with error in units of logioAfQ calculated from multi colour stellar population fitting techniques ; (col. 6) A2800 
Dust correction and error in magnitudes, determined from UV slope fitting ; (col. 7) Effective radius and error in units of kpc from 
Sersic r 1 /" 2D models fits of the Hiqq band data using GALFIT. ; (col. 8) Sersic index and error from Sersic r 1 /" 2D models fits of the 
-H16O band data using GALFIT. ; (col. 9) e-folding star formation time in years calculated from multi-colour stellar population fitting 
techniques (see §2.3). 



Inner Star Formation Growth (IG) : This category is 
for galaxies in which the initial stellar mass added via 
star formation is greater over both regions than the stel- 
lar mass density present, pobserved,inner < pSFAnner and 

pobserved. outer ^ PS F, outer 

In Figure [4] we show examples of the three different 
galaxies classes and in Table [2] we list the numbers of each 
class we have in our sample. 



4 RESULTS 

The profiles for the stellar mass already in place at high 
redshift (§3.1) and the stellar mass added via star formation 
(§3.2) are combined to give an evolved modelled stellar mass 
density profile of the galaxy after evolving for ~ 10 Gyr from 
2 =~ 2.5 until the present day. Using the new stellar mass 
density profiles we fit a new Sersic profile of the same form 
as Equation 5 to examine how the stellar mass added via 
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Figure 4. Example of the three star formation growth classifications. a)Non-significant star formation growth (NG) b)Inner star formation 
growth (IG) c) Outer star formation growth (OG). The blue diamonds represent the observed stellar mass density present at high redshift 
based on the Hi6o-band Sersic profile. The red squares represent the stellar mass density added via star formation to z = 0. The black 
circles represent the combined (evolved) profiles of both the stellar mass density present at high redshift and the stellar mass density 
added via star formation to z = 0. The black dotted line is the best fit Sersic profile to the evolved stellar mass density profile. The 
sub-plot shows the change in the stellar mass density profile from the stellar mass density present at high redshift compared to the 
evolved stellar mass density profile. 



Type 


No. of galaxies 


% of sample 


Non-significant SF Growth (NG) 


29 




Inner SF Growth (IG) 


4 


8-»±2:l 


Outer SF Growth (OG) 


12 


26.7tJ 1 3 1 



Table 2. Evolved massive galaxies using the derived tau model of 
SF evolution separated into the three classifications. Insignificant 
Star Formation Growth (NG), Inner Star Formation Growth (IG), 
Outer Star Formation Growth (OG). We see that nearly half of 
the sample resides in the NG class with a significant fraction in 
the OG class. 



star formation would change the structure and sizes of our 
massive galaxies over time. 



4.1 Stellar Mass 

Figure [3] shows the growth in total stellar mass for all of the 
galaxies within our sample. As stated before the average 
growth for the sample is 91 ± 22%. The evolved total stellar 
masses of our sample of galaxies does not exceed constraints 
placed upon the observed total stellar mass evolution from 
other studies (e.g. Conselice et al. 2007; Brammer et al. 2011; 
Mortlock et al. 2011). This represents the maximal stellar 
mass increase, negating the effects of supernova and other 
types of feedback that would impede star formation and re- 
duce the total amount of stellar mass created. From this 
figure we can also see that there is a clear divide between 
the three classes of galaxies in our sample (see §3.3). The 
NG class galaxies have the smallest change in total stellar 
mass of 22 ± 34%, and do not come close to doubling in 
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stellar mass. The IG class has the largest change in stellar 
mass of 474 ± 89%, and lie exclusively in the top region of 
the figure. The OG class of galaxies in this sample have an 
intermediate mass change of 129 ± 90%. This is a clear seg- 
regation in stellar mass build up via star formation between 
the three classes based upon the star formation locations, 
showing that the three different classes also have differing 
specific star formation rates, with IG galaxies having the 
highest and NG having the lowest. This divide is also present 
in all of the other models of SF we applied to this sample. 



4.2 Structure and Size Evolution 

We find from the Sersic fits to the evolved profiles that the 
average change in n over the whole massive galaxy sample 
is such that the Sersic index goes slightly down, n evo i ve d — 
Tioriginai = An = —0.9 ±0.9. This is consistent with a small 
change with the profile shape over time. 

In Figure [5ji we show that the change in Sersic index 
n differs for the three profile classifications. The NG galax- 
ies lie nearly completely along the 1:1 line, denoting a small 
change from the observed to the evolved n, An = — 0.6±0.1. 
This is expected as these galaxies are classified as having a 
small amount of stellar mass density added via star forma- 
tion compared to the observed stellar mass density. However, 
we note that the NG galaxies with a high original n do not 
fall upon the 1:1, and have a lower n after evolution due to 
small amounts of star formation in the outer regions. 

Also, we find that the majority of the OG and IG galax- 
ies lie below the 1:1 line. This reveals that these galaxies have 
a lowered n value after star formation evolution. Over the 
whole OG class there is a change of An = —1.6 ± 0.4, and 
An = —1.1 ± 0.3 for the IG class. The result is expected for 
the OG class as these are defined as galaxies where there is a 
disparity in the amounts of star formation between the two 
regions, inner and outer. This disparity results in the outer 
regions of the galaxy increasing in stellar mass density, while 
the inner region does not. 

The small changes in Sersic n after star formation evo- 
lution shows that the star formation density within these 
galaxies largely follows the underlying stellar mass density 
profile. In the case of IG galaxies where the stellar mass 
added via star formation dominates over the entire galaxy, 
the initial stellar mass density profile is almost completely 
negligible after evolution but the new profile retains the 
same general shape. 

The effective radius, R e , for our entire sample increases 
by 16 ±5% averaged over the entire sample after our simula- 
tion. Separating the galaxies into our different classifications 
we find that the NG class has a very minor increase in size 
of 4 ± 3%. We find the OG class has an increase in R e of 
37 ± 12%. This increase in the effective radius is due to the 
addition of stellar mass in the outer regions of these galax- 
ies. The IG class has an increase in R e of 36 ± 16%. This 
small increase is most likely related to the non-changing n 
we find for this class. 

Figure [5p shows the evolution of the effective radius be- 
fore and after star formation evolution. We find that galaxies 
in the NG class all lie close to their original effective radius 
with the other two classes having a larger change. We also 
see that the massive galaxies with smaller original effective 



radii have a larger growth in size after star formation evolu- 
tion compared to systems with larger original effective radii. 

4.3 Stellar Migration 

We show above that the star formation within the mas- 
sive galaxies is not sufficient to produced a large growth 
in effective radius. We now investigate stellar migration as 
a method that may also be at work. 

Recent theoretical work suggests that it may be com- 
mon for stars to migrate radically across significant distances 
within spiral galaxies (Sellwood & Binny 2002, Roskar et al. 
2008). These works showed that stellar migration happens 
via processes in the spiral arms of disk galaxies. However, we 
cannot reliably distinguish disk-like galaxies in out sample 
using a Sersic index cut at n = 2.5 because we cannot rule 
out that some of the galaxies with n > 2.5 do not have spi- 
ral like features (e.g. Buitrago et al. 2011). Therefore we add 
into our evolution models a simple stellar migration model 
to all the galaxies in the sample in order to gauge the effect 
this would have on the size and structural changes of the 
evolved galaxy profiles. 

In order to simulate this effect on the galaxies in this 
sample we apply a Gaussian distribution function centred 
on each individual radial bin across each galaxy in the sam- 
ple. This distributes the total stellar mass (the stellar mass 
added via star formation and the in situ stellar mass) across 
the galaxy according to the summed Gaussian distribution: 

M*, mig (R,t)= V M„(i,t)— l -=e~ {J ^- (8) 

where M*(i,t) is the total stellar mass at radius i and at 
time t from equation 6 and a is the width of the Gaussian 
distribution. The Rmin and R ma x are the total range of the 
galaxy radial distribution. Simulations from Roskar et al. 
(2008) showed that the radial migration of stars in a Milky 
Way type disk can change by several kpc over the lifetime of 
the galaxy. To simulate this with some scatter we run a series 
of different widths to represent a wide range of migrations 
with, O.Olkpc < a < l.Okpc. The new stellar mass distri- 
bution is then fit with a new Sersic profile. Figure |S] shows 
the effect of increasing levels of stellar migration upon one 
galaxy at different levels. 

As expected, we find that larger levels of migration have 
an increasing stronger effect on the effective radius. The 
maximal effect on the radius of the galaxy is for the largest 
width Gaussian we applied, l.Okpc. This level of migration 
is of the same order as the effective radii for the galaxies in 
our sample. Therefore to test maximal size growth we use a 
l.Okpc migration from this point on. It is quite possible that 
stellar migration may happen on a larger scale and we exper- 
imented with wider Gaussian distributions but found that 
we could no longer accurately simulate the galaxy evolution 
due to losing stellar mass outside the confines of the simula- 
tion. From the l.Okpc Gaussian we find the effective radius 
grows by 54 ± 19% of the original effective radius. This rep- 
resents an evolved effective radius ~ 3.4 times larger than 
with just the star formation evolution alone. The effect of 
this migration on the Sersic indices on average are consis- 
tent with a small change, with An = —1.1 ± 1.3, similar 
to the non migration case. This is similar to the insignifi- 
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Figure 5. The evolution of the Sersic index (a), and the effective radius (b) via star formation to z = 0. The black circles denote galaxies 
classified as NG. The red triangles denote galaxies classified as OG and the blue squares denote galaxies classified as IG (see §3.3). The 
dashed line in both cases shows a 1 : 1 relation. The typical error bars are shown in the top right hand corner. 



cant change in n we found in the profiles without the stellar 
migration. 

In Figure[7]we show how the Sersic index changes using 
different stellar migration models. When this is applied to 
our different galaxy classes we find that the NG galaxies 
have a increased effective radius of 48 ± 7% over the initial 
effective radius. 

The IG galaxies show a increase in the effective radius 
of 55 ± 15%. Compared to the increased radii from star for- 
mation alone this is a ~ 1.5 times larger result. OG galaxies 
have an effective radius the largest increase with migration 
of 71 ± 18%. This increase in effective radius is ~ 1.9 times 
larger than the non-migration case for the OG class of galax- 
ies in this sample. This is likely due to these galaxies pro- 
ducing more stellar mass in the outer regions than the other 
classes buy definition, and therefore having a larger amount 
of stellar mass already at large radii to move during migra- 
tion. 

In Figure [5] we show how the effective radius changes 
due to the addition of the star formation and stellar mi- 
gration. The galaxy classes that have the highest star for- 
mation rates are affected the most by stellar migration due 
to having more new stellar mass to migrate, with the non- 
changing, non-significant star formation galaxies lying close 
to the non-changing line and the outer and inner star for- 
mation growth galaxies lying above. 



5 DISCUSSION 
5.1 Size Evolution 

As stated previously, recent studies over the last few years 
have found evidence for a dramatic size evolution of massive 
galaxies over the past 10 billion years (e.g Daddi et al. 2005; 
Truijllo et al. 2007; van Dokkum et al. 2010; Buitrago et 
al. 2008, 2011). Current estimates for this size growth argue 
that massive galaxies may grow in size on average up to 
a factor of 3 for disk-like galaxies, while for spheroid-like 



objects this evolution reaches even a factor of 5 since 2 = 3 
(Buitrago et al. 2008). 

In this paper we have shown that the effective radius 
of massive galaxies is altered by the star formation present, 
growing on average by 16 ± 5% from z = 3 to z — 0. This 
value is only ~ 3 — 5% of the total increase in the size of 
massive galaxies from observational studies (e.g. Buitrago et 
al 2008). This indicates that the star formation has a very 
minor contribution to the observable overall size evolution 
at z < 3. 

When we apply a simple model of stellar migration 
to the new stellar mass created via star formation to the 
present day we find that the size of these massive galaxies is 
influenced to a greater extent. The effective radius increases 
by 54 ± 19%. This increase would represent 11 — 18% of the 
total size evolution that massive galaxies undergo between 
2 > 1 and 0. This result shows that the effects of stellar 
mass added via star formation, and any subsequent stellar 
migration, plays a minor role in massive galaxy size evolu- 
tion and only contributes roughly a tenth of the total size 
growth needed to explain the observed size evolution. This 
implies that other evolution mechanisms must also be at 
work to produce the remaining ~ 80% of the observed size 
growth over cosmic time. From also examining the total size 
growth in the other models of evolution (see §5.4.1) we also 
find that the maximal size increase we can obtain can only 
produce ~ 54% of the total observed size growth. 

Recent studies have found that minor and major merg- 
ers have a large influence on the size evolution of massive 
galaxies. These mergers could explain the majority of the 
remaining ~ 80% of the observed size growth unaccounted 
for by the SF via increasing the total stellar mass of the 
galaxies (Bluck et al. 2011). Our results are consistent with 
this view that something other than SF produces the change 
in the sizes of massive galaxies. 



© 2012 RAS, MNRAS 00Q.rnfT8l 



12 Ownsworth et al. 



o 
2 



Q 



10' 2 

10" 

10'° 

10 9 

10 8 

100 
10 

1 



Binned Sersic initiol moss density profile 
Density odded via star formation evolution to z=0 
Moss density profile ofter evoution to z=0 
Fitted evolved Sersic profile 




Q 



R (Kpc) 

(a) 



10' 
10' 
10" 



10 9 

10 8 

100 
10 

1 



Binned Sersic initial moss density profile 
Density odded via star formation evolution to z=0 
Moss density profile ofter evoution to z=0 
Fitted evolved Sersic profile 




10 



R (Kpc) 

(b) 



o 

Cl 



o 



10' 2 

10" 

10'° 

10 9 

10 s 

100 
10 

1 



Binned Migrated Sersic initial mass density profile 
Density odded vio stor formotion evolution to z=0 
Mass density profile after evoution to z=0 
Fitted evolved Sersic profile 



(III. 



10 



R (Kpc) 

(c) 



Figure 6. Example of the effect of the stellar migration models on one example galaxy density profiles with increasing Gaussian widths 
of a = (a)O.fkpc, (b) 0.5kpc, (c) l.Okpc. The non-migration profile can be seen in Figure 4 (c) which is the galaxy we use in this example. 
The sub-plot shows the change in the stellar mass density profile from the stellar mass density present at high redshift compared to the 
evolved stellar mass density profile. 



5.2 Structural Properties 

Recent studies have shown that the massive galaxy popula- 
tion at z ^ 1.5 is dominated by disk like galaxy morphologies 
with n < 2 (e.g., Buitrago et al. 2011; Weinzirl et al. 2011). 
This is in contrast to the local universe where the massive 
galaxy population is almost entirely dominated by spheroids 
(e.g. Baldry et al. 2004; Conselice et al. 2006). This trans- 
formation is also seen through changes in the Sersic index 
of these galaxies from a low value of n at z > 1 to a high 
value of n at z < 1. 

In this study we show that due to the star formation 
present within the massive galaxies at z > 1.5 the Sersic 
index has an insignificant change over cosmic time, An = 
— 0.9T0.9. When we introduce the effects of stellar migration 
to the mass added via star formation the change in Sersic 
index is again negligible over cosmic time with, An = — 1.1± 
1.3. In the other methods of SF evolution we find that the 
change in n is very similar. This implies that with both 



star formation and stellar migration the change to the Sersic 
index is minimal. Also, this does not agree with observations 
of the general increase of n over time. Therefore SF alone 
cannot account for the observed morphological change which 
appear to show that n is increasing over time (e.g Buitrago 
et al. 2011). 



5.3 Spatial Location of Star Formation 

In this study we find that the structural properties of our 
massive galaxies remain largely unchanged after evolution 
via star formation. This unchanging n shows that the loca- 
tion and magnitude of star formation within massive galax- 
ies largely follows the observed initial stellar mass density 
profile. This is most pronounced in the case of the inner 
growth (IG) galaxies. In this class of galaxies the observed 
stellar mass profile is much smaller than the stellar mass 
profile added via star formation. Therefore, for this class 
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Figure 7. The evolution of the Sersic index with due to various stellar migration models, the Gaussian form for stellar migration is 
shown in Equation 8. The values of <r used are (a) O.lkpc, (b) 0.5kpc (c) l.Okpc. The black circles denote galaxies classified as NG. The 
red triangles denote galaxies classified as OG and the blue squares denote galaxies classified as IG. The typical error bars are shown in 
the top right hand corner. 



of galaxy to retain its original Sersic index the stellar mass 
produced via star formation over evolution to the present 
day would have to be produced in amounts which largely 
reflect the already present stellar density i.e. high density 
regions would have a higher star formation rates than lower 
density regions. This was also seen in other ways in Trujillo 
et al.(2007), Buitrago et al. (2008) and Cassata et al. (2010, 
2011). 

The measured ~ 16% growth of the effective radii of 
our massive galaxies due to star formation alone, without 
any stellar migration, reveals that there is star formation 
located in the outer regions of our massive galaxies. This is 
most pronounced in the OG galaxies by definition. In these 
galaxies the surface stellar mass density of the inner region 
remains roughly constant over star formation evolution with 
the outer regions increasing in stellar mass density. Thus in 
our simulated star formation evolution the observed high 
redshift galaxy would become surrounded by an envelope of 



new stellar material over time. With the addition of stellar 
migration this effect becomes more pronounced with newly 
created stellar mass migrating outwards. Recent work exam- 
ining the stellar mass density profiles of high redshift, z > 2, 
and low redshift, z — 0, massive galaxies has shown that the 
density in the core region of low redshift galaxies is compa- 
rable to the density of the compact high redshift galaxies 
(Hopkins et al. 2009; van Dokkum et al. 2010; Carrasco et 
al. 2010). The compact high redshift galaxies have become 
surrounded by an envelope of lower density material from 
z > 2 to 0. This is similar to what we find in the OG class 
of galaxies. 

The models that we use in this study do not account for 
any new gas that can be accreted at later times, at z < 1.5, 
and at early times at z > 3 where we also do not observe 
our sample. This new gas and possible new star formation is 
likely to have a different radial distribution from the current 
in situ gas, with most of the new gas being at larger radii 
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Figure 8. The evolution of the effective radius due to various stellar migration models, the Gaussian form for stellar migration is 
shown in Equation 8. The values of a used are (a) O.lkpc, (b) 0.5kpc (c) l.Okpc. The black circles denote galaxies classified as having 
non-significant star formation growth (NG). The red triangles denote galaxies classified as having outer star formation growth (OG) and 
the blue squares denote galaxies classified as having inner star formation growth (IG). The typical error bars are shown in the top right 
hand corner. 



(Keres et al. 2005; Dekel et al. 2009). Therefore the distri- 
bution of star formation that we observe at high redshift is 
mostly likely the result of previous events of gas accretion 
(see Conselice et al. 2012). However, not all the gas accreted 
may convert into stars immediately, and this gas may re- 
main in the outer portions of these galaxies and may form 
into stars at an epoch later than our observations at z < 1.5, 
which in principle may increase the sizes of these systems at 
a later time, or alter their Sersic indices. 

5.4 Model Limitations 

In this study we have taken a snapshot of our massive galaxy 
sample over 2 Gyr in time, and derived the resulting evo- 
lution based on a derived star formation model. Thus we 
do not take into account any post-observation star forma- 
tion events in our basic model. However this is likely a fair 
assumption due to observations of the majority of massive 



galaxies at z < 1.4 having old stellar populations and red 
colours (e.g. McCarthy et al. 2004; Daddi et al. 2005; Saracco 
et al. 2005; Bundy et al. 2006; Labbe et al. 2006; Conselice et 
al 2007; Mortlock et al. 2011; Griitzbauch et al. 2011). This 
would imply that the SF we observe at z > 1.5 is the last ma- 
jor burst of SF in massive galaxies. The effect of new star 
formation events would increase the total amount of stel- 
lar mass added to the host galaxy. The galaxy's structural 
properties and size could also be affected by these events, 
depending on the location and magnitude of this star for- 
mation as discussed in the previous section. 

Conversely, we also do not take into account any feed- 
back mechanisms that would negatively affect star forma- 
tion rates. Examples of such processes are AGN and super- 
novae feedback. Massive galaxies can spend up to 1/3 of 
their lifetimes in an AGN phase (Hickox et al. 2009, Bluck 
et al. 2011). This phase introduces energy into the interstel- 
lar gas and can expel it from the host galaxy (Schawinski 
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et al. 2006), or heat it such that it cannot cool. Also ongo- 
ing star formation results in the creation of many high mass 
stars which can lose mass during evolution and subsequently 
die in supernovae, thereby lowering the total stellar mass of 
the galaxy. When many supernovae are present in a short 
time the created shock waves introduce vast amounts of en- 
ergy into interstellar gas. The gas can then can be heated or 
ejected from the host galaxy (e.g. Bertone et al. 2007). The 
result of these feedback mechanisms would be a reduction 
of the star formation rate, and the total stellar mass within 
the galaxy would be lower. This decreased amount of stel- 
lar mass added via SF would also result in the stellar mass 
added via star formation having a decreased effect on the 
total size growth and morphological change. 

We also use a very simple models to describe the stellar 
migration that is limited to the extent of the z$*,o band pro- 
files. This means that we can not accurately measure how 
large values of stellar migration would affect the sizes and 
structural properties of our massive galaxies. However even 
though we cannot accurately measure the Sersic index or the 
effective radius of the simulated galaxies with larger values of 
the stellar migration, we find that the stellar mass begins to 
be distributed evenly over all radii, with increasing amounts 
of stellar mass lost outside the confines of the simulation. 
The amount of stellar mass added via star formation moved 
by migration is constant for each galaxy but is distributed 
over wider areas for larger values of stellar migration. This 
results in the stellar mass density added via star formation 
to individual regions of the massive galaxies dropping to 
increasingly smaller values. This implies that with larger 
values of stellar migration, the stellar mass density added 
via star formation would have an increasingly smaller effect 
on the total stellar mass density profile. Therefore even if 
larger values of stellar migration could be simulated in this 
study the change in Sersic index and effective radius after 
star formation evolution and migration would be negligible. 

Stellar migration has also been found, in simulations, 
to be most affected by spiral arms in galaxies (Roskar et 
al. 2011). 73% of the sample of massive galaxies have a 
low Sersic index, n < 2.5, implying a disk-like morphology. 
Within these galaxies we may assume therefore that stel- 
lar migration via disk features may take place, but this is 
far from certain. A few of the galaxies in our sample have 
a high Sersic index, n > 2.5, implying an early-type mor- 
phology, and within these galaxies stellar migration is less 
understood. This does not imply that stellar migration does 
not take place in these galaxies but it must occur by other 
processes than those involving disks. Also, as stated in §4.3 
we cannot reliably distinguish disk-like galaxies in our sam- 
ple using a Sersic index cut because we cannot rule out that 
some of the galaxies with n > 2.5 do not have spiral like 
features (e.g. Buitrago et al. 2011: Mortlock 2012, private 
communication) . 

5.4-1 Evolutionary models 

In this paper we extrapolate the star formation evolution 
using a exponentially declining star formation model based 
on SED derived r values. This value can be uncertain so we 
explore different models of evolution that the star formation 
could follow down to z = 0. Firstly we do not investigate an 
exponentially increasing SFR evolution model because pre- 



vious studies (e.g. Papovich et al. 2011) show that galaxies 
at z < 3 are not well described by this SF history. Therefore 
we investigate the SF evolution models of: constant SFR to 
z — 0, constant SFR to z = 1.5, maximum valid tau and 
minimum valid tau. 



• Constant SFRo to z = 0: This model of evolution as- 
sumes that the massive galaxies we observe at z > 1.5 have 
a very large reservoir of gas and can continue the observed 
SFR over the next lOGyr. This evolutionary method pro- 
duces galaxies in the local universe with very high star for- 
mation rates compared to the galaxies we observe (e.g. Con- 
selice et al 2007). This combined with the fact that over 
the course of their evolution these galaxies will have accu- 
mulated significant amounts of stellar mass with the aver- 
age massive galaxy in this sample increasing its total stellar 
mass by ~ 1500%. This large amount of stellar mass added 
to the galaxies increases the value of R e by 80 ± 20%. This 
is an increase of a factor of 5 over the derived tau model 
in effective radius growth, but still only 16 — 27% of the 
observed size evolution. This model of evolution is highly 
unlikely due to the many features of this model that we do 
not observe in the local universe, such as very large stellar 
mass growth leading to very massive galaxies with stellar 
masses over 10 13 Mq (e.g. Brammer et al. 2011; Conselice et 
al 2011; Mortlock et al. 2011, all find that the stellar mass 
growth at the massive end of the luminosity function is on 
the order of 200% from z > 1.5 to 0) and very high star 
formation rates of 100's of solar masses per year. 

• Constant SFRo to z — 1.5: This model of SF evolu- 
tion is based on the observation that the majority of mas- 
sive galaxies at z < 1.4 have old stellar populations and 
red colours (e.g. Conselice et al 2007: Mortlock et al. 2011, 
Griitzbauch et al. 2011). This would imply that these galax- 
ies have turned off their SF before z = 1.5. To model this we 
employed a constant observed SFR until z = 1.5 at which 
point the SFR is reduced to 0. In this evolution scenario 
the total stellar mass of the massive galaxies is increased by 
126 ±20%. The effective radii in this model are increased on 
average by 37± 19%. This is a factor of ~ 2.5 larger then the 
increase from the derived tau model. This is still insignifi- 
cant to the total observed size increase. This model has a 
very similar effect on the change in n, An — — 1.1 ± 1.1, as 
the derived tau model. 

• Maximum valid tau to z = 0: In this model of evolution 
we use the largest value of tau derived for our galaxy sample, 
r = 2.71 x 10 9 yr. We apply this exponentially declining rate 
to all the galaxies in the sample. In this scenario we obtain a 
large average increase in total stellar mass of the sample of 
377 ± 172%. The change in the effective radii of this model 
is on average R e — 57 ± 33% a factor of ~ 3.8 larger than 
the derived tau model of evolution. This increase in effective 
radius is still only ~ 11— 19% of the observed size evolution. 
The change in n for this model, An = —1.5 ± 1.7 is similar 
to change for the derived tau model. 

• Minimum valid tau to z = 0: This model is similar to 
the previous model. Except that the minimum valid tau, 
r = 1.2 x 10 8 yr, is used to extrapolate the SF. This would 
give the shortest time scale that the SF would occur. In 
this model the average galaxy in the sample increases its 
stellar mass by only ~ 17%. This very small increase in mass 
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is accompanied by an equally small change in R e , average 
AR e = 3 ± 1%, and n, An = -0.4 ± 0.6. 

From this investigation of different models of SF evolu- 
tion to z — we find that the value in the increase of the 
effective radii of the massive galaxies can at no point fully 
explain the total observed size increase. The valid models of 
SF evolution that we applied can only produce a factor of 
~ 3.8 times larger than the size increased we obtained from 
using the derived tau model at maximum. The change in 
Sersic index in all the models are within the error consistent 
with the answer obtained from the derived tau model used 
in this paper. 



5.4-2 Dust Gradients 

In this paper we assume that the dust obscuration is con- 
stant across the radius of individual galaxies. From studies 
of local and distant studies this may not be the case. Colour 
gradients in the local universe have been shown to corre- 
spond to age and dust gradients (e.g. Boquien et al. 2011; 
Smith et al. 2012). 

We apply a dust gradient to our sample of massive 
galaxies that allows the attenuation due to dust to vary 
within the given error across each galaxy. This is done in 
two ways. A positive dust gradient with higher attenuation 
towards the outer regions of the galaxy, and a negative dust 
gradient with higher dust attenuation towards the central 
regions of the galaxy. 

In the positive gradient case we find that the average 
increase in the effective radius was 68 ± 36% larger than the 
original measured effective radius. This is a factor of ~ 4.5 
larger change than the growth in R e we obtain from using 
a radially constant dust correction. From this gradient the 
change in n is largely the same as before but with a much 
larger scatter, An = —0.9 ±2.0. The positive gradient could 
contribute a maximum of ~ 23% to the 300 — 500% size 
growth. 

In the negative gradient case we find that the average 
increase in R e is minimal, Ai? e = 7 ± 3%. This small in- 
crease in the effective radius is accompanied by a change in 
n that is very similar to most other cases, An = —1.0 ± 1.0. 
This negative gradient case would seem to produce a very 
small increase in the effective radii of our sample and only 
contribute a maximum of ~ 2% to the total observed size 
growth. 

Neither of the gradient cases that we applied to the 
sample are able to fully explain the observed size growth or 
observed change in Sersic index. 



within each galaxy. We also construct new stellar mass pro- 
files of the in situ stellar mass at high redshift from Sersic 
fits to rest-frame optical, i/i6o-band, data. We combine the 
two stellar mass profiles to produce an evolved stellar mass 
profile. 

We then fit a new Sersic profile to the evolved profile, 
from which we examine what effect the resulting stellar mass 
distribution added via star formation has on the structure 
and size of each individual galaxy. In summary: 

• We find three different profiles of star formation within 
the massive galaxies in this sample, Non-significant Star For- 
mation Growth (NG), Outer Star Formation Growth (OG) 
and Inner Star Formation Growth (IG) (see §3.3). With most 
of this sample of massive galaxies falling in to NG class using 
the derived tau model of evolution. 

• We find that the star formation we observe at high red- 
shift, and its effects on galaxy sizes, is not large enough to 
fully explain the observed size evolution of effective radius of 
~ 300 — 500%. Star formation alone alone can only produce 
an increase in effective radius on the order of ~ 16% over 
the whole sample. This value can vary as much a a factor 
of 4.5 by using different evolution mechanisms but is always 
insufficient to fully explain the observations. 

• We find that over the whole sample of massive galax- 
ies the stellar mass added via star formation has a slight 
effect on the Sersic index of the evolved galaxy profile such 
that they decrease. This indicates that the star formation 
within these galaxies follows the same radial distribution as 
the original stellar mass profile. This also implies that star 
formation evolution has a minimal effect on structural evo- 
lution between z ~ 3 and the present day. 

• The increase in effective radius can be enhanced by 
adding in the effects of stellar migration to the stellar mass 
created via star formation. This increases the total effective 
radius growth to ~ 55%, which is still however much smaller 
than the total observed size increase. 

We conclude that due to the lack of sufficient size 
growth and Sersic evolution by star formation and stellar 
migration other mechanisms must contribute a large propor- 
tion to account for the observed structural evolution from 
2 > 1 to the present day. Recent studies by Bluck et al. 
(2011) have found that minor and major mergers have a 
large influence on the size of massive galaxies possibly con- 
tributing the remaining 80% of size growth needed to explain 
the observed trends. Large surveys such as CANDELS and 
future telescopes such as JWST and E-ELT will provide the 
quality of data that is required to explore the star forma- 
tion locations of lower mass galaxies and probe resolved star 
formation at higher redshifts for similarly massive galaxies. 



6 SUMMARY 

We investigate the resolved star formation properties of a 
sample of 45 massive galaxies (M* > 10 1 Mq) within a red- 
shift range of 1.5 ^ z ^ 3 detected in the GOODS NICMOS 
Survey, a HST J/i6o-band imaging survey. We derive the 
star formation rate as a function of radius using rest frame 
UV data from deep zsso ACS imaging. The star formation 
present at high redshift is then extrapolated to z — 0, and 
we examine the stellar mass produced in individual regions 
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